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ABSTRACT 

We present the first detection of the mass donor star in the prototype X-ray binary Scorpius X-f . 
Phase resolved spectroscopy revealed narrow emission lines components from the irradiated secondary 
star. Radial velocity fits to the Bowen blend emission are used to establish an absolute orbital ephemeris 
and derive an accurate value for the systemic velocity; 7 = —113.8 ± 0.6 km s -1 . Emission from the 
irradiated front side of the secondary leads to solid limits to the radial velocity of the mass donor of 
87 < K<2 < 148 km s _1 . In conjunction with an upper limit to the velocity of the primary K\ < 53 km 
s -1 , we derive a firm limit on the mass ratio of Sco X-l of q < 0.61. A likely set of system parameters 
satisfying the various constraints as well as the recent inclination estimate by Fomalont et al. (2001) is 
K\ = 40 km s _1 , K2 — 133 km s _1 , q — 0.30 corresponding to component masses of A/i = 1.4M , 
M 2 = O.42M for an orbital inclination of 38°. 

Subject headings: accretion, accretion discs — stars, individual (Scorpius X-l) — X-rays: binaries 



1. INTRODUCTION 

Scorpius X-l is the prototype low-mass X-ray binary 
(LMXB), consisting of a neutron star accreting matter 
from a Roche lobe filling companion star. It is the bright- 
est persistent X-ray source in the sky and has been the 
target of detailed analysis for the last 35 years in all en- 
ergy bands, from radio to X-rays. An orbital period of 
18.9 hr was reported by Gottlieb, Wright & Liller (1975) 
from brightness eye estimates on 1068 plates during the in- 
terval 1890-1974. The optical B-band light curve shows a 
0.13 mag amplitude modulation (Augusteijn et al. 1998) 
which is being interpreted as the visibility of the X-ray 
heated companion. Radial velocity curves of Hell A4686 
and HI confirmed the orbital period and indicated that the 
inferior conjunction of the emission line regions are close 
to the photometric minimum and, hence, they must origi- 
nate near the compact object (e.g. LaSala & Thorstensen 
1985). Bradshaw et al. (1999) measured the trigonomet- 
ric parallax of Sco X-l using VLB A radio observations and 
deduce a distance of 2.8 ± 0.3 kpc. Intensive radio obser- 
vations were also used to measure the properties of the 
twin radio lobes around the source which are inclined at 
an angle of 44 ± 6° relative to the line of sight (Fomalont 
et al. 2001). Despite intensive efforts, spectroscopic signa- 
tures of the companion star have not yet been observed in 
optical-IR data (e.g. Schachter, Filippenko & Kahn 1989; 
Bandyopadhyay et al. 1997,1999). Here we present high 
resolution spectroscopy of Sco X-l which enable us to de- 
tect the companion star for the first time. 

2. OBSERVATIONS AND REDUCTION 

The phase resolved spectroscopy of Scorpius X-l was 
obtained with the 4.2m William Hershel Telescope on La 
Palma in combination with the ISIS double beam spectro- 
graph. The R1200B grating centred at 4601A delivered a 



wavelength coverage of 4150-5050A at 0.45A/pixel on the 
2048x4096 pixel EEV CCD detector. A 1.2" slit width 
resulted in a spectral resolution of 0.84A, corresponding 
to 55 km s _1 at our central wavelength. The data were 
part of a campaign to look for circular polarisation from 
accreting binaries, and to this end a quarter- wave plate 
and a Savart calcite slab were placed in the optical path 
(see Tinbergen & Rutten (1992) for more details). The slit 
position angle was adjusted to remain close to parallactic 
throughout the observations. A total of 137 spectra were 
obtained in the course of three consecutive nights covering 
75% of the 18.9 hour orbital cycle. A journal of observa- 
tions can be found in Table 1. 

The insertion of the calcite slab in the light path results 
in the projection of two target beams on the detector. The 
intensity ratio then represents the fraction of circular po- 
larisation. For the analysis presented in this Letter, we 
ignored the polarisation characteristics of the light, and 
combined the extracted spectra from both beams. Frames 
were first de-biased using a median bias frame obtained 
by combining 10-15 individual bias exposures. The over- 
scan strip was then used to subtract any residual bias from 
each frame. A median of 29 Tungsten exposures was con- 
structed at each night in order to perform flat field cor- 
rection. Each beam was then optimally extracted. Some 
frames were affected by transient CCD pickup noise pro- 
ducing variable bias structure at the level of 3-5 ADU. 
For those frames affected (around 8 on the second night 
and 10 on the third night), we introduced an additional 
reduction step before optimal extraction that subtracted 
the pickup noise using a 2nd order polynomial fit to each 
CCD row. Regular CuAr arclamp exposures allowed us 
to establish an accurate wavelength scale for each spec- 
trum. Finally, a wide slit exposure of the flux standard 
BD28+4211 (Oke 1990) was used to perform flux calibra- 
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tion using the MOLLY data analysis package. The two 
extracted spectra for each frame were then added to pro- 
duce the final set of 137 calibrated Sco X-l spectra. 

3. ANALYSIS 

3.1. The emission line spectrum 

The optical spectrum of Sco X-l is dominated by strong 
and broad emission lines from the Balmer series as well 
as Hel/Hell. However, closer inspection revealed a large 
number of weaker and very narrow emission line features 
that showed significant Dopplcr motions as a function of 
orbital phase. In particular, the strong Bowen blend near 
4640A was resolved into an underlying broad component 
as well as several narrow components from the individ- 
ual Bowen transitions. Broad and highly variable emis- 
sion around 4640A was already noted by Sandage et al. 
(1966) in the paper discussing the identification of the op- 
tical counterpart of Sco X-l. A more detailed investiga- 
tion of the same data by Westphal, Sandage & Kristian 
(1968) did also indicate complex and strongly time de- 
pendent fine structure in the broad emission structures 
and some preliminary evidence for velocity changes in the 
strong Balmer and Hell lines. More recently, Bowen emis- 
sion in Sco X-l was noted by Willis et al. (1980) while 
Schachter, Filippenko & Kahn (1989) discuss the Bowen 
fluorescence process in Sco X-l in detail. However, nei- 
ther had the spectral resolution to resolve the individual 
components. We will see in the next sections that these 
narrow emission features can be attributed to the irra- 
diated companion star, thus revealing a direct signature 
from the mass donor in Sco X-l for the first time. In Fig- 
ure 1 we present the average normalised spectrum of Sco 
X-l in the rest frame of the companion star to facilitate 
the identification of the narrow features. Most lines were 
identified using the planetary nebulae emission line lists 
from Kaler (1976) as well as the Kurucz atomic line list. 
Many narrow lines from Oil arc present as well as a few 
Silll, CII and NIII transitions. 

3.2. Bowen emission from the donor 

In order to establish the origin of the narrow emission 
line components, we measured the velocity of the narrow 
Bowen components as a function of phase. As can be seen 
in Figure 1, the strongest Bowen components are the two 
NIII components at 4634.13A and 4640.64A, and the CIII 
components at 4647.4/4650. lA. In Figure 2, we display 
the Bowen blend and HeII4686 emission in the form of a 
trailed spectrogram. All the individual spectra were first 
normalised to the continuum using a 3rd order polynomial 
fit. The two NIII components as well as the CIII compo- 
nent at 4647A clearly move in phase with each other, and 
are extremely narrow. In the background more diffuse 
emission can be seen to move roughly in anti-phase with 
these sharp components. For comparison, the broad and 
weakly double peaked Hell line shows completely different 
line kinematics. 

The velocities of the various components were measured 
using multiple Gaussian fits to the individual spectra. At 
first we fitted the Bowen profiles using 3 narrow Gaus- 
sians corresponding to the 3 strongest Bowen components 
as well as an underlying broad component. The FWHM 
of the Gaussians were determined from a fit to the average 



Bowen blend profile and fixed to 80 km s _1 for the narrow 
components and 1550 km s _1 for the broad component. 
The position and strength for each Gaussian was a free pa- 
rameter for each individual spectrum. The derived radial 
velocities of the sharp components were identical within 
the error bars. We thus decided to reduce the number 
of free parameters in our fits and give the three narrow 
Gaussians a common radial velocity, but an independent 
strength. The 6 parameters were optimised for each ob- 
served spectrum through least square fitting to determine 
the common velocity of the three narrow components and 
its associated formal error. 

Figure 3 plots the derived radial velocity curve of the 
Bowen components together with a least squares sinusoidal 
fit. The best fit gives a velocity amplitude of 77.2 ±0.4 km 
s _1 with a mean velocity of — 113.8±0.6 km s _1 . Given the 
extreme narrowness of the components, they cannot orig- 
inate from the accretion flow or hot-spot that is expected 
to be responsible for the broad emission line components. 
Such narrow emission components are, however, regularly 
seen in the line profiles of cataclysmic variables where they 
originate from the irradiated front side of the companion 
star (e.g. Harlaftis et al. 1999, Schwope et al. 2000). This 
has also been observed in the intermediate mass X-ray bi- 
nary Her X-l (Still et al. 1997), where narrow Bowen blend 
components from the mass donor are present. An origin on 
the mass donor is furthermore supported by the fact that 
the crossing of the emission components from blue shift to 
red shift occurs near minimum light, and that the strength 
of the emission is weakest then, since we are looking at the 
unirradiated backside of the secondary. The inferred sys- 
temic velocity of 7 = —113.8 ± 0.6 km s _1 furthermore 
agrees well with previous estimates using the strong emis- 
sion lines (e.g. Lasala & Thorstensen 1985). We thus have 
established a radial velocity curve of the mass donor in Sco 
X-l for the first time. The radial velocity curve directly 
gives us an unbiased value for the systemic velocity of the 
binary and establishes the correct absolute orbital phases. 
From the sinusoidal fit we established the following orbital 
ephemeris which we will use throughout this paper; 

To(HJD) = 2451358.568(3) + 0.787313(1)^ 

where we took the value for the orbital period as de- 
rived by Gottlieb, Wright & Lillcr (1975). Using our 
ephemeris, the phase of minimum light as established by 
Gottlieb, Wright & Liller corresponds to an absolute phase 
of (j) = 0.04±0.04, where we have allowed for the uncertain- 
ties in both To and P. This confirms that the photometric 
modulation is indeed related to the X-ray heated hemi- 
sphere of the donor star. To values based on the inferior 
conjunction of the Hell A4686 emission, on the other hand, 
are delayed by 0.10±0.02 (Cowley & Crampton 1975) and 
0.08 ± 0.02 (LaSala & Thorstensen 1985), indicating that 
the emission line regions trail the compact object. We will 
revisit these emission line velocities in Section 3.3. 

Since we have detected narrow lines from the mass donor 
in emission rather than absorption, they must originate 
from the front half of the Roche lobe. The derived radial 
velocity amplitude is thus biased towards lower values de- 
pending on the detailed distribution of the line emission 
over the Roche lobe surface. This in analogy with the bias 
towards larger velocities when using absorption lines that 
are predominantly formed near the backside of the Roche 
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lobe (e.g. Wade & Home 1988). The fact that the ra- 
dial velocity curve is not strictly sinusoidal also suggests 
that the emission is not distributed uniformly. However, 
allowing for an ellipticity in the fit does not significantly 
improve the x 2 of the fit. We can thus constrain the true 
velocity of the centre of mass of the mass donor (-K2) to be 
K2 > 77.2 km s _1 . Without any knowledge of the mass 
ratio and the extent of disc shadowing effects in Sco X-l, 
it is very difficult to quantify the K-correction that needs 
to be applied to the Bowen blend velocities in order to de- 
fine the true K2 . A more detailed modelling of the phase 
dependence of the emission from the mass donor, both in 
terms of velocity as well as strength is required, which will 
be discussed in a future paper. 

3.3. The radial velocity of the primary 

Direct inspection of the trailed spectra (Figure 2) clearly 
shows that the Hell A4686 line is anti-phased with the 
CIII/NIII narrow emissions and, therefore, it must come 
from near the compact object. Significant phase shifts be- 
tween the line core and centroid have been reported, due 
to possible contamination by complex variable hot-spots 
(e.g. Cowley & Crampton 1975). Therefore, we have fol- 
lowed the double-Gaussian technique to estimate the ve- 
locity semi-amplitude of the compact object K\ from the 
line wings (see Schneider & Young 1980 for details). We 
employed a two Gaussian band-pass with FWHM=100 km 
s _1 and Gaussian separations a = 200 — 1000 km s _1 in 
steps of 100 km s _1 . Radial velocity curves of different line 
sections yield consistent results with velocity semi ampli- 
tudes in the range 50-56 km s _1 and blue-to-red crossing 
phases 0.64-0.60. Phase delays of 0.1 — 0.2 relative to the 
expected motion of the compact object, are typically ob- 
served in cataclysmic variables (CVs), probably due to 
residual hot-spot contamination. We have also attempted 
to extrapolate the velocity points towards phase 0.5 using 
the light-centre technique (see Marsh 1988 for details) but 
the method fails given the very small derived phase drift 
as we move from the line core into the wings. Applying 
the same technique to the H/3 line produces similar values 
for the line core, however the phase offset increases dra- 
matically in the line wings. We therefore choose to use 
the Hell lines as our basis for an estimation of K\. The 
bottom panel of Figure 3 presents the radial velocity curve 
of Hell A4686 for a = 600 km s _1 , before velocity points 
start to be corrupted by continuum noise. A formal sine- 
wave fit provides K = 53 ± 1 km s _1 which we take as an 
upper limit to the true K\. 

3.4. Doppler mapping 

Many of the narrow lines that can be identified in Fig- 
ure 1, move in phase with the secondary star and thus 
are also the result of reprocessed emission from the irradi- 
ated front side of the Roche lobe. Unfortunately, due to 
weakness of the lines, a Gaussian fit to each spectrum is 
not possible. Instead, we use Doppler tomography (Marsh 
& Home 1988) to measure the radial velocity amplitudes 
of the many weak narrow emission lines that are present. 
By mapping the observed data onto a velocity coordinate 
frame, Doppler mapping makes use of all data at once, 
and can thus be used for features that are too weak to 
be separated in each individual spectrum. One effectively 



resolves the distribution of line emission in the corotating 
frame of the binary system, providing an excellent tool for 
identifying the origin and kinematics of the various emis- 
sion components. It also allows us to map the emission 
line distribution from the broad HcI,HcII and Balmcr lines. 
Secondary star emission is readily identified in Doppler to- 
mograms since the solid body rotation of the Roche lobe 
is mapped to a corresponding Roche lobe area along the 
positive V y axis. We employed a maximum entropy im- 
plementation of Doppler tomography whereby the data is 
fitted under the added constraint of maximising the en- 
tropy, i.e. image smoothness of the tomogram (Marsh & 
Home 1988). This reduces the presence of noise artifacts 
in the recovered tomograms and allows a simultaneous fit 
to a number of heavily blended lines. 

To check that the method of measuring the velocity of 
the narrow emission lines using tomography is consistent 
with the method of employing (multiple) Gaussian fits, we 
started by producing tomograms for the Bowen blend com- 
ponents. The Bowen blend data was fitted simultaneously 
with three Doppler images corresponding to the rest wave- 
lengths of the three strongest Bowen emission components. 
The tomograms revealed the secondary star emission as a 
sharp spot along the V y axis, and a 2D Gaussian fit was 
used to accurately measure the position of the spot. This 
lead to inferred velocities of 76 ± 1 , 76 ± 1 and 77 ± 1 km 
s _1 respectively, confirming that this method can be used 
to extract reliable radial velocities. We then repeated the 
same procedure for a range of other weak emission lines 
present in the spectrum of Sco X-l. The inferred velocities 
for emission from the mass donor spanned from 66 km s _1 
for Silll at 4568 A, 78 km s" 1 for Oil doublet at 4415 A to 
87 km s _1 for NIII at 4379A. An example tomogram show- 
ing secondary star emission from the NIII line at 4514A is 
plotted in Figure 4. The fitted velocity based on this to- 
mogram gives 79 km s _1 . Unfortunately, the width of the 
narrow components are of the order of our instrumental 
resolution and so are not well resolved. However, since 
only the front half of the Roche lobe can contribute to the 
line emission, we can confidently increase our lower limit 
to K2 > 87 km s _1 . This technique in conjunction with 
higher spectral resolution data could be used to map the 
emission on the Roche lobe as a function of ionisation po- 
tential such as done in for example Harlaftis (1999), and 
provide even more stringent constraints to the true value 
of K 2 . 

In Figure 4 we also present Doppler tomograms of the 
strong and broad Hf3, HeII4686 and HeI4921 lines for com- 
parison. Clearly their line dynamics are radically different 
and are dominated by contributions from the accretion 
disc around the primary. However, significant asymme- 
tries are present in all lines, leading to localised hot-spots 
in the tomograms. In the case of Hell at 4686A, we can 
identify a clear emission contribution from the secondary 
star, whereas the Balmer and Hel lines do not show a 
secondary star emission component. If anything, the H/3 
tomogram indicates a depression of line flux at the loca- 
tion of the mass donor. The tomograms directly illus- 
trate that the use of emission line velocities as a measure 
of the radial velocity amplitude of the compact object is 
severely affected and distorted by asymmetries in the line 
profiles. Indeed, our analysis in Section 3.3 revealed signif- 
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icant phase offsets between the phasing of the secondary 
star components and the broad emission lines. The differ- 
ent phase offsets that exist between the velocities of the 
H/3 line a versus Hell are due to the fact that the Hell is 
mainly distorted by emission from the mass donor, whereas 
the H/3 line is distorted by a spot in the lower left quad- 
rant of the Doppler image. The fact that our data reveals 
the same phase shifts as found by Lasala & Thorstensen 
(1985), suggest this asymmetry is a persistent feature in 
Sco X-l. 

The complex morphology of the Doppler tomograms 
complicates the reliable estimation of K\ . If the emitting 
gas was symmetrically distributed around the primary, the 
Doppler image would reveal a circular distribution centred 
on V x — 0, V y = —K\. In order to estimate the true cen- 
tre of the emission line distribution, we subtracted the 
symmetrical component from the H/3 Doppler maps and 
inspected the residuals. We cycled the assumed centre of 
symmetry between -100 and +25 km s _1 in steps of 2 km 
s _1 and measured the mean of the residuals at high ve- 
locities in the upper right quadrant which is not affected 
by the spot. The mean residuals show a pronounced mini- 



mum around an assumed centre of symmetry of V y 



-40 



km s 1 . This suggests that the true value of K\ is indeed 
closer to 40 km s _1 for Sco X-l. 

4. DISCUSSION 

We have presented the first direct detection of the mass 
donor in Sco X-l by resolving emission components orig- 
inating from the irradiated front side of its Roche lobe. 
Radial velocity fits were used to establish an absolute or- 
bital ephemcris and measure the systemic velocity to high 
accuracy; 7 = —113.8 ± 0.6 km s _1 . A lower limit to 
the radial velocity of the mass donor is given by the ra- 
dial velocity amplitude of the sharp Bowen blend emission 
components which leads to K2 > 77 km s _1 . Detection of 
secondary star emission at even higher velocities in other 
transitions lead to a more stringent lower limit of K^ > 87 
km a~ 1 . The strong broad emission line show orbital mo- 
tions roughly in anti-phase with the secondary star. We 
confirm significant phase shifts between the radial veloc- 
ity curves of the strong lines and the absolute ephemeris. 
Doppler tomography reveals that these phase shifts are 
due to persistent asymmetries in the emission line distri- 
bution around the compact object. The formal double 
Gaussian method leads to a firm upper limit to the pri- 
mary radial velocity amplitude of K\ < 53 km s _1 . A 
more reliable value is derived by searching for the optimal 
centre of symmetry in the accretion disc emission which 
leads to an estimated K\ value of 40 ± 5 km s _1 . Our 
absolute ephemcris confirms previous suggestions that the 
phase of optical minimum light corresponds closely to or- 
bital phase zero, when the contribution from the heated 
mass donor is smallest. 

The combination of our upper limit to K\ and at the 
same time a lower limit to K2 leads to a firm limit to 
the mass ratio of q = M2/M1 < 0.69 when using the for- 
mal values of K\ < 53 and K^ > 77 as derived from the 
Gaussian fits. Using the more restrictive values K\ ~ 40, 
K2 > 87, we derive a limit of q < 0.46. Another con- 
straint is given by the fact that the observed emission 
originates from the front half of the Roche lobe and its 



velocity thus lies in between that of the true K2 value 
and the velocity of the LI point (Vli). Since we observe 
emission at velocities as low as 66 km s" 1 , we can thus 
derive an upper limit on K2 for a given mass ratio in or- 
der to make sure that the velocity of LI is not greater 
than 66 km s _1 . In Figure 5 we plot the various limits 
on K2 as a function of mass ratio. As can be seen in the 
plot, the combination of Vli < 66 and K2 < K\/q leads 
to K2 < 148 km s" 1 regardless of the mass ratio. This 
upper limit proves interesting from the point of view of 
the recent inclination estimate based on the active radio 
lobes of Sco X-l (Fomalont et al. 2001). If the orienta- 
tion of the radio lobes corresponds to the rotation axis of 
the orbital plane, we can estimate the component masses 
using Kepler's law. With a typical neutron star mass of 
Mi ~ 1.4M©, this inclination (44 ±6°) predicts relatively 
high K2 values across the allowed mass ratio range (Figure 
5). This suggests that a slightly lower orbital inclination 
seems more likely unless the neutron star is atypical and 
significantly less massive than the Chandresekhar limit. 
However, given the uncertainty on the inclination, all con- 
straints are still compatible with a 1.4 Mq neutron star 
within lcr. A likely set of parameters for Sco X-l given the 
various constraints could be K\ = 40 km s _1 , K2 = 133 
km s _1 , q = 0.30 corresponding to component masses of 
Mi = 1.4M Q , M 2 = 0.42M Q for an orbital inclination of 
38°. Given the orbital period of Sco X-l, the companion 
star must then be a significantly evolved subgiant having 
a radius of more than twice of the corresponding main 
sequence radius of ~ O.5i?0. Its mean density is solely de- 
termined by the orbital period which gives p — 0.31g cm~ 3 
compared to ~ 5.6 g cm~ 3 for a 0.42M Q main sequence 
star. We note that these parameters are also consistent 
with Sco X-l being a persistent X-ray binary, with an ap- 
proximate mass transfer rate of M2 ~8x 10~ 10 MQyr~ 1 
using the prescription from King, Kolb & Burberi (1996). 
This is enough to maintain the irradiated-disc hotter than 
6500 K everywhere avoiding global instabilities and thus 
outbursts. 

Bandyopadhyay et al. (1997,1999) fail to detect any 
features from the companion star in their IR spectra of 
Sco X-l at two epochs. Most notably, the CO band-head 
near 2.2/xm that is present in late type stars is absent. 
However, the appearance of the donor star will be affected 
by irradiation and both epochs were taken near orbital 
phase 0.5 (using our derived ephemeris) when these effects 
would be particulary severe. In addition, the use of the 
CO-bands for spectral type classification can be mislead- 
ing because their strength is strongly dependent on effec- 
tive temperature and, furthermore, the C,0 abundances 
can be significantly non-solar for companion stars in in- 
teractive binaries (e.g. Howell 2001). In order to derive 
more accurate system parameters, we need to be able to 
model the K-correction, and use higher spectral resolution 
data to resolve the emission across the Roche lobe for a 
range of transitions using the methods described in this 
paper. Detecting the companion star in absorption should 
be possible around true orbital phase zero when the effects 
of X-ray heating are minimal. 
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Table 1 
Journal of observations 




Date 


UT interval 


Number of spectra Exposure time 


Orbital phase interval 


28/6/1999 
29/6/1999 
30/6/1999 


21:55-02:10 
20:55-02:15 
20:56-01:55 


40 300s 
51 300s 
46 300s 


0.81-1.03 
0.03-0.31 
0.30-0.56 
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Fig. 1. — The average normalised spectrum of Sco X-l in the rest frame of the donor star. The three panels are on the same intensity scale 
with the strong HcII4686 and H/3 going off-scale in order to highlight the weaker lines. The peak strength in Hell is around 1.46 compared 
to the continuum, whereas it is 1.39 for H/3. 
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Fig. 2. — Trailed spectrogram of the Bowcn blend and HcII4686. Three sharp components moving in phase with each other stand out from 
the broad Bowen blend emission. Small gaps separate the data from the three nights. 
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Fig. 3. — Top: radial velocity curve of the narrow Bowen blend emission components as derived from a multiple Gaussian fit. Data is 
repeated over two cycles. Best sinusoidal fit is over plotted as a solid line. The orbital phase of the binary is calculated using the best fit 
radial velocity curve to determine the zero point. Bottom: Radial velocity curve of the Hell emission line based on a double Gaussian fit with 
separation 600 km s _1 and FWHM of 100 km s _1 . Sinusoidal fit is over plotted indicating a phase shift of 0.6 with respect to the donor star 
emission. 
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Fig. 4. — Doppler tomograms reconstructed from the H/3, HcI4921, HcII4686 and NIII 4514 lines using maximum entropy optimisation. 
Both grayscale as well as contours span the full scale of the maps. The origin corresponds to the centre of mass of the binary system using 
7 = -113.8 km s _1 . 
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Fig. 5. — Plotted are constraints to the value of K2 as a function of mass ratio. The horizontal solid line is the lower limit to K2 as derived 
from the observed line emission (K2 > 87). The dot-dashed lines denote the upper limits on K2 given that K2 < K\/q for K\ = 53 (upper 
curve) and Ki=40 km s~ 1 . The requirement that the velocity of the LI point may not exceed 66 km s~ 1 leads to the dashed curve. The 
two dotted curves denote the values of K2 if the primary is a 1.4 Mq neutron star. The two curves correspond to ±1<t limits for the orbital 
inclination based on the orientation of the radio lobes (44 ±6°). 



